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Abstract 

The classical limit for generalized partition functions is obtained using coher- 
ent states. In this framework it is presented a general procedure to obtain all 
the corrections to the classical limit. In particular, the first and second order 
quantum corrections are worked out explicitly, and the classical limit for the 
Tsallis thermostatistics is determined. The results of this work generalize the 
ones obtained by E. Wigner (Phys. Rev. 40 (1932) 749) for usual statistical 
mechanics. 

PACS number(s): 05.20.-y, 05.70.-a, 05.30.Ch 



Typeset using REVTgX 



*Permanent address: Departamento de Fisica - Universidade Estadual de Maringa 
Av. Colombo 5790, 87020-900 Maringa, PR, Brazil 

1 



I. INTRODUCTION 



Long-range interactions are present in many physical contexts. The nature of these in- 
teractions can be connected to spatial and temporal dependencies. For instance, anomalous 
diffusion astrophysics with long-range (gravitational) interactions some magnetic 
systems [PHHT2], some surface tension questions In such cases we have a nonextensive 

behavior. This fact indicates that usual thermodynamics and statistical mechanics deserves 
some changes. Thus, it is important to generalize the concepts based in the usual thermo- 
statistics. In this direction, it has been recently analyzed the Legendre transform structure 
15| , [T6[| and stability conditions |16|] in a very general context. The present work is devoted 



for kind of study. More precisely, the purpose of this work is to obtain the classical limit 
and its quantum corrections for arbitrary partition functions. 

It is important to obtain the classical limit and its corrections because, in some cases, it is 
easier to calculate the classical partition function than the quantum one. This problem was 



solved first for the usual statistical mechanics by Wigner |17j using the Wigner functions, and 
by Kirkwood employing the Bloch equations (see also reference |T9| , |20[| ). The procedure 
employed by Kirkwood can not be applied easily for generalized partition functions. On 
the other hand, the procedure used by Wigner can be in principle extended for arbitrary 
partition functions, but in this work we prefer to employ coherent states |?TJ. By using 
coherent states we classify and give the prescription to calculate all the possible corrections 
for the classical partition function. Employing this prescription we obtain explicitly the first 
two corrections. 

The general results obtained in this work can be in principle applied to an arbitrary 
thermostatistics. In particular, they are used in the context of Tsallis thermostatiscs. This 
thermostatistics was proposed in order to cover nonextensive systems (long-range micro- 
scopic memory, long-range forces, fractal space time). Thus, the Tsallis thermostatistics 



| 22| has interesting proprieties and has been applied in many situations, like for instance, 
self-gravitating systems p3|-[25|1, two-dimensional-like turbulence Levy-like p6|-p0|1 and 
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correlated-like [|3T|-|3"4}1 anomalous diffusion, solar neutrino problem 
ory [^6| and magnetic systems |37|-[4(J. 
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, linear response the- 



II. QUANTUM CORRECTIONS 

In this work we consider a partition function for a very general statistical mechanics, 

Z = Trf(H), (1) 

where Tr represents the trace, H is the Hamiltonian of the system, and f(x) is the gener- 
alization of the exponential function of the usual thermostatistics. 

To calculate the classical limit of ([!]) and its corrections it is convenient to employ 
coherent states. In particular, two basic properties of the coherent states are important for 
our calculations. The first one is the definition of coherent states, namely 

a n | z >= z n | z > , (2) 

where 



m n uj i \ 

and \ z >= Yin=i I z n >■ I n this expression for | z >, iV is the dimension of the config- 
uration space, and each | z n > is normalized, < z n \ z n >= 1. Moreover, q n and p n are 
respectively the classical values for coordinate and momentum n. The other property is the 
overcompleteness relation, 

f f-f dp n dq n I I /c .x 

1= yn-2^ri*><*i- ( 5 ) 

Furthermore, employing the previous notation, we use A to represent an operator and A to 
represent its classical analog. 



To obtain the corrections to the classical partition function we basically expand f(H) 
around its classical value, f{H). More precisely, we write H(p n , q n ) as H(p n + 5p n , q n + 5q n ) 
and express f(H) as a power series in 5H, 

f(H + 5H) = f(H) + §j*H+ - + - . , (6) 

where 

SH = H(p n + 5p n , q n + 5q n ) - H(p n , q n ) , (7) 

N -2 

H = H(p n ,q n ) = Y,^ + V(q n ), (8) 

n=l Am n 

8q n = q n — q n and Sp n = p n — p n . Note that in the present analysis we are considering a 
Hamiltonian that depends only on the spatial coordinates and its corresponding momenta. 
By using the relation (|5|) we can now express the partition function as 

z = ffl*g±<*\f(H + m\'>, 0) 

n=l 

and to calculate it we employ the expansion presented above. The first term of this expansion 
gives the classical partition function, 



'0 



N 



dp n dq n 



n^f/w- do) 



The corrections to ( |T0D can be classified using a simple but important result, i.e. 

< z | bX x ...bX u | z > oc h u/2 , (11) 

where 5X n can be 5q n or 5p n . In fact, the expression fllTD comes directly from the definition 
of 5q n and 5p n , 



^ = Vsb( <5a; ' +<te ") (12) 



and 



fipn = i J mn ^ (Sal - 5a,nj , (13) 



where 5a r 



and Sal = al — z*. Furthermore, when u is an odd number a direct 



n 



calculation leads to < z \ 5Xi...5X u \ z >= 0. Therefore, the corrections to ( |I0D are 
proportional to integer powers of h. Thus, the partition function can be written as 



z = z + nz 1 + n 2 z 2 + ... . 



(14) 



For instance, the first and second quantum corrections, hZi and HZ 2 , comes respectively 
from the second and fourth order expansions in 5q k and 5pk- 

Let us now proceed to calculate the first quantum correction. First we note, by means 
of equation (§, (@) and flU), that 



< z < 



' 5q n 5q k ' 




' 1/ (m n uj) ' 


Sqjpk 


> z >= h < 


i 




> 


Spjq k 




—i 


, SpnSpk . 







5nk 



(15) 



Moreover, from (|6]), ([j]) and (§) these mean values are present only in 5H and (5H) 2 . 
Consequently, by using the equations (0),(||) and (|1^) we conclude that 

1 



< z I 5H I z > 



2 ^ 

Z n=l 



+ 



d 2 V 



2 2m„w 2 



+ 0(h 2 ) 



(16) 



and 



<z\{5Hf\z> = ^-Y, 

Z n=l 



dV 



It follows directly from (| 



m„c<j 2 \<9g r , 
dH), (|16D and (0) that 



C(/i 2 



(17) 



d# TTljdH 2 



d 2 V df (dV\ d 2 f 
7)qJdH + [dq~ dlP 



;is) 



j=l I V"-- 1 - 1 ""J u '- l - t / ""r 

This expression can be written in a more compact form. In fact, after some derivations we 
verify that 



7 _ w f TT dqndp 

n=l j=l 



E 



1/ + j d _( d _L^_ 

8pj \ (III J </•/_, V </. /( 



dH ) rrijU! 2 dqj \dqj dH 
Furthermore, the expression (|T9|) is more convenient than (]T8|) for future discussions. 
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(19) 



The calculation of the second quantum correction, h 2 Z 2 , is similar to the previous one. 



In fact, to obtain the analogous formula to ( fLq) 



< 5X n 5X k 5X s 5X u > = -m; a " /2 mr t/2 m; as/2 m7" /2 (w)" ( '" +fft+a '" 1/2 



(20) 



it is necessary only the equations (0), ( |T2| ) and (pT3|). In the expression ( p0|) it was used a 
compact notation suggested in the Eq. (JfTJ) because there are sixteen combinations of 5qj 
and Spj. In this notation 5Xj = h l ^ 2 i^ l ~ a ^/ 2 {rrijUj)^' T ^ 2 {5ajj + o\,<5a) is equal to 5q.j when 
ex,- = 1 and equal to 5pj when <jj = —I. By using the expressions (fL5|) and (|20f) we can 
obtain after some calculation the contributions proportional to ft and h 2 in < z \ {H) n \ z >, 



i. e. 



< z iff z >= — > 
11 2 ^ 

n=l 



1 i a 2 v 

2 2m n uj 2 dq 2 



h 2 1 r) 4 V 

Au^Sm n m k dqldy 1 J 



(21) 



h, , N 

< z\(5H) 2 \z>=^Y: 

A n=l 



ft 2 N 

i,fc=i 



4 I uj 2 mjm k 



dV d 3 V 
dqj dqjdql 



lcFVd 2 V 1 / <9 2 V y 

4 <9g| 0g£ 2 \dqjdq k ) 



5j k 



+ 



d 2 V 



' 2m.,/ dg 2 . 



+ ^Q + ^)KO(ft 8 ), (22) 



< z | («f) 3 |z>4E 



4 ,1^-, I uj 2 rrijm k 



1 /<9^\ 0V dFdV d 2 V 
+ 



PjPfc <9 2 V 1 fdV 
3rrijm k dqjdq k 2m j \dq 3 



2\dq j J dq\ dqjdq k dqjdq k 

2 

i_ 
2m 



+ 



|fc- = l+A fc + 2)[ + 0(fi 3 



p 2 d 2 y 

2mjm k dq\ 
(23) 



and 



J j^ =1 \4:UJ 2 m j m k \dqj J \dq k ) 2m j m k \dq k ) itu , itu k 



+ o(h 3 ) 



(24) 



It is convenient now to group the terms of the previous expressions with the same power 
in uj. Thus, the expression for Z 2 can be written as 
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where 



tt dq n dp Tl 



1 d 4 V \ df 



+ 



1 N 

- y 

8 mjirik y\A dq 2 dq\ J dH 

1 (dV\ 2 d 2 V dVdV d 2 V 
+ 



2\dq j ) dqi dqjdqkdqjdqk 



;h + h 



w 2 L 



dV d 3 V ld 2 Vd 2 V 



dqj dqjdql 4 dq] dq\ 
d 3 f 



1 

2 



' d 2 V 
dqjdq k 



dH 3 



+ 



1 (dV\ 2 (dV\ 2 



4 \dq 3 



dq k/ 



dH 4 



(25) 



dH 2 



(26) 



and 



1 N 

8 j,k=i 



. rrii 



2m i 



+ 2m j (<%) ' 6kk 



d 2 V 
dH 3 



dH 2 



+ 



d 2 v PjPk d 2 v 



2m^m k dq\ 3mjm k dqjdq k 



+ 



p 2 



2m j m k \dq k ) 
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N 

E 

j,k=i 



'1 SpA (f£_ 
4 + T/ dlP + 



^ kk + 2) 



dH 4 



dH 3 \2mj2mk) dH 4 ) 



(27) 



(28) 



As in the Z 1 case it is more convenient, for the following discussions, to write Ii, I 2 and I 3 
in a more compact form, 



h 



N 

E 



Q3 



1 dV df 



32 dqjdql \mjm k dqj dH t 



(29) 



and 



4 fri m fe ^ dq 2 ) dH 2 6 ,^ 

N / ™ , \ 2 



AT 



Pi Pfc d 2 V 1 



ayy d 4 f 1 " 



mjm k dqjdq k mjm k \dq k J 
' p k dV , Pj rf 2 / 



dH 3 



m k dqj 2m j dq k ) dH 2 



N 



d 



^3 qoEa, 



32 ^9pj L Pj '^ 2 



N 8 ( d 2 f\ 



(30) 



(31) 



This general procedure can be extended to other orders in H to obtains further quantum 
corrections, but we prefer to discuss the previous results. 

In Z\ and Z 2 are present terms which depends of the frequency u employed in the 
definition of the coherent states. Note that these terms are surface terms. However, a 
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classical limit can not depend on the particular choice of the coherent states basis. This 
apparent contradiction can be eliminated if we consider that in a macroscopic body the 
surface effects can be neglected when we compare them with the bulk contributions. By 
using this consideration we conclude that the hZ\ contribution to the partition function can 
be neglected, and only I 2 is relevant in the computation of h 2 Z 2 . Therefore, when these 
considerations are employed we verify directly that 



N 

n 

n=l 
1 N 

- y 

j,k=l 



dp n dq n 



f(H)+h 2 j n 



dp n dq n I 

2irh ' ' " J 2vr^ | 

Pj Vk d 2 V + _J_ fdvV 
mjm k dqjdq k mjTn k \dq k ) 



A' 



1 (d 2 V\ d 2 f 



4 ^ m k \ dq 2 k J dH 2 



d 3 f 



N 



dH^ S^ymkdqu) dH* 



PkdvY d*f 



0{h 2 ) . (32) 



Note that, as it is expected, f(H) = exp(—j3H) in the previous expression recovers the usual 



results 17.19.2 



III. THE CASE OF TSALLIS THERMOSTATISTICS 



We apply now the previous findings for Tsallis thermostatistics p2| . The entropy is 
defined as 

S, = ^-f"" , (33) 

where q is a real parameter and p is the density matrix. Note that q gives, essentially, the 
measure of the non-extensivity of (^). Furthermore, to thermal averages, one has to use 
the g-expectation value defined as 

< A > q = Tr(piA) , (34) 

for some relevant observable A. Thus, the canonical distribution is given by 

f(H)={l-(l-q)PH}^ . (35) 

When (1 — q)(3E n > 1, where E n are the eigenvalue of H, the probability is not be positive 
defined quantity. To avoid this behavior it is usually assumed that f(H). In particular, the 
Tsallis thermostatistics reduces to the usual one in the limit q — ► 1. 
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The substitution of (|35f ) into the previous results gives the partition function for Tsallis 
thermostatistics up to h 2 order. In this case we have 




nq — (n — 1) 

= (-P) n q(2q - l)(3g - 2) ... [(n - l)q - (n - 2)][1 - (1 - q)/3H] ^ . (36) 



The factor (3 n above allow us to arrive at an important conclusion. In fact, in the limit (3 — > 
the corrections to the classical limit disappear, i. e., the quantum corrections to the classical 
partition function are irrelevant in the high temperature limit. In this case we are supposing 

nq — (n — 1) 

basically that the derivatives [1 — (1 — q)f3H] ~ in (|32|) do not diverge. Note that the 
above conclusion is a direct consequence of the fact that the dependence of the function 
f(H) with j3 and H appears only in the form f3H. Thus, the previous conclusion can be 
applied not only to the Tsallis thermostatistics but to a wide class of possible partition 
functions. 



A. The harmonic oscillator 



To exemplify the application of the above expressions for a concrete case, we consider a 
harmonic oscillator whose Hamiltonian is H = p 2 /(2m) + mu 2 x 2 . For this one-dimensional 
case, (^) is reduced to 
dpdx 



2nh 



^ 9 r dpdx [ 



dpdx f uo 2 d 2 f 1 ( uo 2 



p 2 + mio^x 2 



d 3 f 



to 



2 2 
-X p 



dm] 



+ 0(h 2 ), (37) 



4 dH 2 6 \ m' J dH 3 8 

where the derivatives are obtained through (^6[). In fl3"TD we have to evaluate the integrals 
for two cases: q < 1 and q > 1. 

Let us consider first the case of q < 1. Since the integrands appearing in ([J7|) are 
composed of successive derivatives obtained from (|36|) , the relevant integrals are of the form 



dxdy 1 — (ax 2 + by z 



(38) 



dxdyx 2 



(ax 2 + by" 



na 



-3/2^-1/2 



2(l + a)(2 + <x) 



(39) 



and 

/ dxdyx 2 y 2 \l - (ax 2 + by 2 )Y = 7ra" 3/2 r 3/2 — --^ — -- , (40) 

J L y ! \ 8(1 + a)(2 + a)(3 + a) ' V ; 

with a > — 1 in order to have well defined integrals. The above integrals are more easily 
performed if the transformation x = cT l l 2 r cos 6 and y = 6 -1 / 2 rsin0 is employed. By 
substituting the resulting integrals in (|37D we obtain for the partition function 

One observes that in the limiting case q —>■ 1 we recover the usual result. Moreover, the sec- 
ond factor which accounts for the quantum corrections in the first order in Ti is independent 
of q for the case of harmonic oscillator. It is important to stress that the above results are 
meaningful for q > 2/3 in order to achieve the convergence of all the integrals in fl37|) . In 
the general case, the lower limit for q is fixed by the higher order derivative of f(H) present 
in the expansion for Z. 

The calculations for the case q > 1 are similar to the previous one and lead to the same 



expression for Z, i.e., (|4"I). However, in this case the upper limit is q < 2, and it is now 



fixed by the exponent of the first term in the expansion for Z. 



IV. CONCLUSIONS 

In this paper we have presented for the first time a general procedure to obtain all the 
quantum corrections for generalized partition functions. As underlined before, this procedure 
is a generalization of the Wigner expansion for usual statistical mechanics. However, instead 
to use Wigner representation, we have used the coherent states representation, which is 
convenient for our purposes. The general formalism was used to explicitly evaluate the 
first quantum correction to the classical case. This result was applied to a non-extensive 
(Tsallis) thermostatistics. As an application we have considered the case of the harmonic 
oscillator. In particular, we have shown that the dominant term in the partition function 
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depends on q. It is reduced to the usual one when q — > 1 as expected. Moreover, the first 
quantum correction for the partition function of the harmonic oscillator coming from the 
Tsallis thermostatistics is shown to be the same as the usual one. We have also determined 
the specific heat of the system for q > 1 and q < 1. 
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